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I.  INTRODUCTION  AND  EQUATIONS  OF  MOTION 


A  one  dimensional  continuum  with  one  strain  variable  and  one  internal 
variable  may  be  used  to  model  either  a  long  slender  rod  or  longitudinal 

1  2 

deformation  in  a  granular  material.  *  In  a  rod  the  strain  is  the  usual  axial 
strain,  and  in  a  granular  material  it  is  the  overall  longitudinal  strain.  The 
Internal  variable  measures  radial  strain  in  the  rod  and  introduces  radial 
forces  and  shears  into  the  equations  of  motion.  In  the  granular  material  it 
measures  the  volume  fraction  of  solid  material  and  introduces  forces 
associated  with  the  symmetrical  expansion  or  collapse  of  pores.  The  internal 
variable  also  gives  rise  to  an  extra  Inertial  effect,  which  in  the  rod  is  due 
to  radial  motion  and  in  the  granular  material  is  due  to  pore  dilation. 

The  equations  that  govern  these  two  cases  may  be  set  down  side  by  side  for 
comparison.  Only  a  purely  mechanloal  version  of  either  theory  will  be 
considered . 


Rod 

Granular  Material 

Linear  Momentum 

S'  *  pw 

o*  *  P; 

(1) 

Micromomentum 

Q'  -  P  *  Jjp*2u  , 

h'  ♦  g  s  pkv 

(2) 

Kinematics 

2  »  Z  ♦  w 

v  s  velocity 

(3) 

e  *  w' ,  v  «  w 

e  s  strain 

r  •  R(1+u) 

v  s  volume  fraction 
of  granules 

q  *  u' 

q  *  v' 

Stored  Energy 

W(u,q,e) 

e(v,v» ,e) 

(4) 

Stress  Potentials 

s  «  we 

o  *  eE 

(5) 

h  *  eVf 

g  x  -e 

2T.  W.  Wright ,  " Nonlinear  Waves  in  Rods ,"  in  Proceedings  of  the  IUTAM  Symposium 
on  Finite  Elasticity ,  D.  E.  Carlson  and  R.  T.  Shield,  eds.,  Martinus  Nijhoff 
(1981). - 

2«T,  W.  Runziato  and  S.  C.  Cowin,  "A  Nonlinear  Theory  of  Elastic  Materials  with 
Voids."  Arch.  Rat.  Mech.  Anal.,  72,  17 S  (1979). 


Here  Z  and  t  are  independent  variables  for  space  and  time,  the  prime  and  dot 
denote  partial  differentiation  in  the  usual  way,  and  subscripts  denote  differ¬ 
entiation  with  respect  to  the  arguments  of  the  function.  In  both  cases  the 
stored  energy  is  an  even  function  of  its  second  argument.  This  latter  fact 
has  a  profound  effect  on  the  types  of  waves  that  can  occur . 

Clearly  the  two  continuum  models  are  mathematically  identical  even  though 
the  physical  interpretation  of  the  various  terms  is  dissimilar.  In  the  rod 
the  material  point  with  cylindrical  coordinates  R,Z  is  located  instantaneously 
in  space  at  r,z  so  the  independent  kinematic  variables  are  the  axial  displace¬ 
ment  w  and  the  radial  strain  u,  which  are  assumed  to  depend  only  on  Z,t.  In 
the  granular  material  v  and  e correspond  exactly  to  w  and  w*  ,  and  the  volume 
fraction  v pi ays  the  same  role  that  u  does  for  the  rod.  In  the  rod  p is  the 
reference  density,  and  a  is  the  initial  radius.  In  the  granular  materials  is 
also  the  reference  density,  but  of  the  bulk  material,  not  just  the  granules, 
and  k  is  called  the  equilibrated  inertia,  which  gives  rise  to  an  extra  com¬ 
ponent  of  kinetic  energy  due  to  the  local  effects  of  pore  dilation.  In  the 
rod  S  is  the  axial  engineering  stress,  P  is  the  sum  of  the  average  radial  and 
circumferential  stresses,  and  Q  is  the  average  radial  moment  of  radial  shear 
stress.  In  the  granular  material  a  is  the  longitudinal  stress,  g  is  the 
intrinsic  equilibrated  body  force,  and  h  is  the  equilibrated  stress.  For 
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further  details  the  reader  is  referred  to  the  original  papers  *  but  since  all 
stresses  are  derived  from  similar  potentials,  it  is  evident  that,  except  for 
the  signs  of  P  and  g,  the  analogy  between  the  elastic  rod  and  the  granular 
material  is  complete. 

On  the  basis  of  this  analogy  alone  an  important  fact  emerges.  Since  it 
is  clear  on  physical  grounds  that  the  radius  a  is  the  natural  length  scale  for 
all  motions  of  the  rod,  the  quantity  /2k  is  the  natural  length  scale  for  all 
motions  in  the  granular  material.  The  remainder  of  this  paper  will  use 
terminology  appropriate  for  the  rod  model,  but  naturally  all  results  will 
apply  equally  well  to  the  granular  material. 


II.  ONE  DIMENSIONAL  SHXKS 


The  integral  form  of  a  general  one  dimensional  conservation  law  may  be  written 
as  follows. 

^  (Z,t)  dZ  ■  [£G]  b(t)  +  [m]  b(t^  s  dZ  •  (6) 

(t)  a(t)  a(t)  J  a(t) 

The  end  of  points  a(t)  and  b(t)  are  allowed  to  move  arbitrarily  with  respect 
to  the  material,  f(Z,t)  is  the  conserved  quantity,  [•]  signifies  the  differ¬ 
ence  of  the  bracketed  quantity  between  end  points,  G  is  the  speed  of  the  end 
point  relative  to  the  material,  m  represents  the  flux  of  f  added  at  the  end 
points,  and  s  is  the  Interior  supply  of  f.  In  words  the  equation  simply  says 
that  the  rate  of  ohange  of  some  field  quantity  within  an  arbitrary  interval  is 
equal  to  the  sum  of  convective  transport  through  the  ends  of  the  Interval,  the 


flux  of  f  at  the  ends,  end  the  supply  of  f  throughout  the  interior  of  the 
interval.  In  the  usual  way, it  is  postulated  that  this  equation  holds  even  if 
there  is  a  discontinuity  in  the  field  variables,  say  at  a  moving  point  with 
position  c(t)  in  the  interior  of  the  interval.  Consequently,  as  is  well 
known,  the  Integral  form  of  the  balance  law  yields  both  the  differential  form 
and  the  jimp  conditions  for  shock  waves. 

£  ■  m*  ♦  s  (7) 

V  [£]  ♦  [m]  -  0  (8) 


where  V  =  c(t)  is  the  shock  speed  and  £•]  =  (*)”  -  (*)+  is  the  jump  across 
e(t). 

Since  it  is  reasonable  to  assume  that  the  rod  diameter  cannot  change 
abruptly,  u  must  be  continuous.  Conservation  of  mass  yields  nothing  new,  and 
conservation  of  linear  mcmentun  with  f  =  pv  and  m  s  S  gives  the  familiar  jump 
condition 


pV  [v]  ♦  [S]  =  0 


(9) 


but  conservation  of  micromomentum  with  f  s  1/2  pa  u  and  m  a  Q  yields  an  un¬ 
familiar  Jump  condition 


1/2  pa2  V  [u]  +  [Q]  *  0. 

With  the  required  compatibility  conditions  added 

[v]  *  V  [e]  ■  0  ,  [lij  ♦  V  [qj  »  0 
the  jump  conditions  become. 

PV2  ■  |fj-  and  1/2  pa2  V2  «  • 

3 

These  conditions  were  given  previously  by  Nunziato  and  Walsh.  In  that  same 
paper  they  also  asserted  that  q  (\>x  in  their  terms)  must  be  continuous  in  a 

strong  shock  wave.  In  the  next  section  it  will  be  shown  that,  in  fact: 

1)  There  are  two  types  of  shock  wave  possible; 
il)  Although  q  may  be  continuous  across  one  type,  in  general  it  is  not; 
and 

ill)  In  the  second  type  of  shock  wave  q  is  always  discontinuous,  and 
although  e  is  usually  discontinuous  as  well,  [e]  will  generally  be  weaker  than 
Cq]. 


7 

J.  W.  Nunziato  and  E.  K.  Walsh,  "One-dimensional  Shook  Waves  in  Uniformly 
Distributed  Granular  Materials,"  Int.  J.  Solids  Structures ,  14,  681  (1978). 


(10) 

(ID 

(12) 


The  first  type  of  shock  is  the  familiar  longitudinal  shock  for  both  the 
rod  and  the  granular  material.  The  second  type  is  not  so  familiar,  however. 

In  the  rod  it  is  associated  with  an  abrupt  change  in  u' ,  so  the  surface  of  the 
rod  will  show  an  abrupt  change  of  slope  at  the  shock.  In  the  granular  mater¬ 
ial  there  will  be  an  abrupt  change  in  the  gradient  of  solid  volume  fraction, 
v'.  It  is  legitimate  to  call  this  second  type  of  discontinuity  a  shock  since 
it  involves  jinaps  in  derivatives  of  order  one  lower  than  those  that  appear  in 
the  corresponding  differential  equation. 


III.  WEAK  SHOCK  WAVES 


The  jumps  in  the  gradients  of  w  and  u  may  be  expressed  as  follows. 

e"  ■  e+  ♦  a,  q'  »  q+  +  6  .  (13) 

With  these  expressions  substituted  in  (12)  the  jump  conditions  become 

PV2  a  *  W  (u+,  q+  +  6,  e+  +  a)  -  W  (u+,  q+,  e+)  (14) 

fc  t 

l/2pa2V2S  *  W  (u+,  q+  +  B,  e+  +  ct)  -  W  (u+,  q+,  e+)  .  (IS) 

H  4 


Since  these  two  equations  have  the  three  unknowns  (a,  0,  V),  it  is  to  be 
expected  that  solutions  will  be  one  dimensional  curves  in  a  three  dimensional 
space.  The+exact  shape  and  location  of  these  curves  will  depend  parametri¬ 
cally  on  (u  ,  q+,  e+).  Clearly  one  branch  is  always  given  by  the  trivial  case 
a  =  8  =  0  with  V  arbitrary,  but  bifurcations  may  occur  where  the  Jacobian  of 
equations  (14)  and  (15)  taken  with  respect  to  o  and  0  vanishes  at  a  =  B  s  0. 

J  (o,6)  -  W£e  -  pV2 


J  (0,0)  =  0 

Equation  (17)  determines  the  values  of  V  at  which  bifurcations  occur.  In 
fact,  it  is  just  the  characteristic  condition  for  acceleration  waves  (see 
equations  (3.22)  -  (3-24)  in  Reference  4).  To  look  for  other  branches  near 
the  bifurcation  points,  other  than  the  trivial  branch,  first  expand  the  energy 
W  in  powers  of  a  and  6  in  equations  (14)  and  (15).  In  the  simplest  case 
u+  =  q+  =  e*  =  0  .  Bifurcations  in  (16)  will  appear  at 

pV2  •  W  (0,0,0)  and  l/2pa2V2  =  W  (0,0,0)  (18) 

MM 


J •  W.  Nunziato  and  E.  K.  Walsh ,  "On  the  Influence  of  Void  Compaction  and 
Material  Non-uniformity  on  the  Propagation  of  One -dimensional  Acceleration 
Waves  in  Granular  Materials,"  Arch.  Rat.  Mech.  Anal. .  64,  299  (1977). 


Recall  that  W  is  an  even  function  in  q  so  that  odd  derivatives  in  q  will  not 
appear.  The  expansion  now  appears  as  follows. 


PV2.  .  W%,  .  1/2  W°w  B2  .1/2  p2  .  1/2  «  82  . 

- < 

*  1/6  C/  * . 

(19) 

1/2  p.2  V2  6  -  <Jq8  .  »“qeBa  .  1/6  ^  B3  .  1/2  H^e.2  . . 

•  •  (20) 

The  superscript  0  indicates  that  the  derivatives  of  W  in  (19)  and  (20)  are  to 
be  evaluated  at  u+  =  q+  =  e+  =  q  •  It  is  possible  to  regard  either  a  or  B  as 
the  primary  variable. 

Case  i)  To  examine  the  branch  near  the  first  of  the  bifurcations  in  (18)  let  B 
2 

and  V  be  given  terms  of  a  as  a  power  series. 


V2  »  ♦  VjO  +  VjQ2  +  . 

2  3 

B  •  BjOi  ♦  ®2a  *  63a  +  * 

Substitution  of  (21)  and  (22)  into  (19)  and  (20) 
lowest  powers  of  a  shows  that 


• •  •  •  (21) 
•  •  •  • 

(22) 

and  comparison  of  terms  in 


Bj  ■  0,  i  »  1,2,3, . 

■  <. 
pvl  ■ 

»V2  ’  «« . 


(23) 

(24) 

(25) 

(26) 


This  is  the  case  discussed  by  Nunziato  and  Walsh,  More  generally,  if  |aq^«  l, 
but  +  0  ,  and  u+  ^  0,  e+  /  0  as  well,  terms  in  all  powers  of  a  and  B  will 
appear  throughout  the  expansions  (19)  and  (20).  If  account  is  taken  of  the 
fact  that  odd  derivatives  of  W  with  respect  to  q  are  odd  in  q  ,  the  net  result 
is  that  (23)  -  (26)  will  all  be  modified  by  additional  terms  of  order  0(aq+) 
or  higher.  It  is  easily  worked  out  that  the  lowest  order  corrections  are 

PV0  "  W«  +  °<aV2)  (27) 


W  q 

qqc 


l/2a  W  -W 
cc  qq 


♦  0(a3  q+3) 


(28) 
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In  (27)  and  (28)  tha  superscript  ♦  indicates  that  the  arguments  of  the 
derivatives  of  W  are  (  u  ,  0,  e  ). 

Case  ii)  X?  examine  the  branch  near  the  second  of  the  bifurcations  in  (18) 
let  a  and  V2be  given  in  terms  of  8  as  a  power  series.  Consider  th^  simplest 
case  first  with  q*  =  0.  It  turns  out  that  the  series  for  a  and  V  1  only  have 
even  powers  of  6  in  this  case. 


I 


where  =  W*£  /  p,  c2  =  VQ,as  in  (36)  above.  In  the  rod  the  speeds  c1  and 

Cg  are  associated  with  bulk  longitudinal  and  shear  speeds  [1],  In  the 

granular  material  c1  is  associated  with  the  compressibility  of  the  granules 

without  voids,  and  c 2  with  changes  in  the  void  fraction  without  compression  of 

the  granules. 4  Ordinarily  c  >  c„  ,  but  that  is  not  necessarily  true  for  all 

values  of  u+  and  £+.  The  point  is  that  both  V2  and  V4  may  be  either  positive 

or  negative  depending  on  material  properties  and  the  strain  state  ahead  of  the 
wave. 


As  before,  a  more  general  case  occurs  if  0-*jaq+  |«1.  Then  the 
expansions  in  (29)  and  (30)  will  contain  all  powers  of  8  .  The  new 
coefficients,  that  is  the  coefficients  of  the  odd  powers  of  $  in  the  equations 
that  replace  (29)  and  (30),  will  all  be  of  order  0(aq+)  and  the  ones  already 
calculated  from  (3D  -  (35),  that  is,  the  coefficients  of  the  even  powers  of  B, 

2  +2 

will  change  only  by  order  0(a  q  ).  This  result  again  follows  from  the  fact 
that  W  is  even  in  q. 


u 


u 


IV.  SHOCK  STABILITY 


Shock  speeds  have  been  calculated  above  simply  as  possible  solutions  of 
(1*1)  and  (15).  It  has  been  convenient  to  give  the  speeds  parametrically  in 
powers  of  the  amplitude,  using  either  a  or  B as  appropriate  for  the  type  of 
wave.  However,  it  is  not  reasonable  to  accept  all  values  of  shock  speed 
computed  in  this  manner  without  some  consideration  of  shock  stability.  Since 
it  is  not  possible  to  appeal  to  a  thermodynamic  argument  in  a  purely 
mechanical  theory  such  as  this,  it  is  necessary  to  adopt  some  other  criterion 

such  as  the  Lax  stability  criterion, ^  whereby  it  is  required  that  the 
characteristic  speed  behind  (ahead  of)  the  shock  of  the  same  family  is  greater 
than  (less  than)  the  shock  speed.  Clearly  a  necessary  condition  is  that  the 
shock  speed  for  nonzero  amplitude  must  be  greater  than  the  bifurcation  speed. 
This  condition  is  met  if  the  shock  speed  is  a  monotonically  increasing 
function  of  the  amplitude;  but  without  detailed  calcuation,  it  does  not  seem 
possible  to  say  whether  the  characteristic  speed  behind  the  shock  will  be 
greater  than  the  shock  speed.  Some  of  the  possibilities  that  may  occur  are 
shown  in  Figures  1  through  5. 

Figure  1  shows  the  case  discussed  in  Reference  3.  With  q+  r  0,  the  solid 
line  shows  the  common  example  of  compressive  shocks.  The  dashed  line  is  the 
continuation  of  the  solution  to  (14)  and  (15),  but  now  this  branch  is  not 
allowed  on  stability  grounds  as  discussed  above.  If  q+  ^  0,  the  picture 


5P.  Lax,  Hyperbolic  Systems  of  Conservation  Laws  II, "  Com.  Pure  Appj.  Math 
10,  537  (1957). 
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a)  b) 


Figure  1.  a)  Sketch  of  shock  speed  vs.  amplitude  and  b)  Jump  in  q  vs.  ampli¬ 
tude  for  quasi- longitudinal  shocks.  The  solid  and  interrupted  lines 

show  allowed  shocks  for  q+  »  0  and  q+  j*  0  respectively.  The  dashed 
and  dotted  lines  show  forbidden  shocks  for  the  same  cases. 


14 


changes  only  slightly  with  the  interrupted  line  showing  allowed  shocks  and  the 
dotted  line  showing  forbidden  shocks.  Figure  1b  shows  that  8  =  0  in  the 
simplest  case,  but  that  in  general  g  £  0  when  q+  f  o. 

The  possibilities  for  the  second  type  of  shock  are  much  richer.  Some  of 
these  are  shown  in  Figures  2-5  where  it  has  been  assumed  throughout  that 

c1  f  c2«  Figure  2  shows  the  case  when  V2  <0  and  v2  >0.  If  q+  =0.  no  shocks 

at  all  are  possible  unless  the  amplitude  is  greater  than  a  critical  value,  and 

then  the  amplitude  may  be  of  either  sign.  If  q+  t  0  ,  then  shocks  of  one  sign 
only  are  possible  for  a  limited  range  of  amplitude,  but  for  large  enough 
amplitudes  shocks  of  either  sign  are  possible.  The  fact  that  q*  f  0  has 
merely  introduced  some  asymmetry  into  the  curve.  Likewise  the  curve  fora  as  a 
function  of  B is  symmetric  if  q+  =  0  and  asymmetric  if  q+  f  o  .  Figures  3-5 
show  the  other  possibilities.  Note  that  in  every  case  stable  shocks  may  exist 
for  some  range  of  amplitude  if  q+  j*  0,  and  in  nearly  every  case  if  q+  =  0. 


V.  STEADY  WAVES 


Steady  waves  in  a  nonlinear  elastic  rod  have  been  discussed  previously 
for  the  special  case  when  the  stored  energy  function  can  be  decomposed  into 
W(u,q,e  )  =  (u.o  )  +  V.^ (q)  •  Here  the  analysis  for  steady  waves  will  he 

briefly  reviewed,  but  for  an  arbitrary  stored  energy.  In  a  steady  wave  all 
field  variables  depend  only  on  a  special  combination  of  Z  and  t,  namely 
5  *  Z  -  ct  where  c  is  a  constant  speed  of  propagation.  The  partial 
differential  equations  (1)  and  (2)  become  ordinary  differential  equations,  and 
it  was  shown  in  Reference  1  that  the  solution  can  be  reduced  to  quadratures  of 
the  following  form, 

e  -  f(u;A,B,c)  (39) 

u'2  -  g(u:A,  B,c)  (40) 


where  the  functions  f  and  g  satisfy  two  integrals  of  the  motion  identically. 
These  are 

F (u,g,f)  =  Wf  -  pc2f  -  A  =  0  (41) 

G(u,g,f)  =  2W  g  +  Af  -  W  -  l/4pc2a2g  +  l/2pc2f2  -  B  =  0  (42) 

o 

where  A  and  B  are  constants  of  integration.  A  power  series  solution  for  f  and 
g  in  terms  of  u  may  be  found  by  differentiating  (41)  and  (42)  with  respect  to 
u  and  solving  each  successive  pair  of  equations  for  the  derivatives  of  f  and 
g.  The  choice  of  initial  values  for  f,  g,  and  u  determines  the  constants  A 

and  B.  For  example,  choosing  f  (0)  =  g  (0)  =  0  leads  to  values  of  A,B,  and 

the  derivatives  at  u  =  0  as  follows. 
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v2-  V 


a) 


Figure  4.  a)  Sketch  of  s 
tude  for  T&di&l 


1/4pc-.-)  ♦  2  Wuff  £u  ♦  Wfff  r 

♦**_  (*...  ♦  w..,  f) 


l/4pc  a")  (W 


f  -  pc  )  (W  -  1/4  pc  a  ) 


(4Vu  +  3V  *uu  *  Wuu« 


♦2W  -f  +  W  ,,f  *  W  „f 


(Wg  -  1/4  pc  a  ) 


Truncated  versions  of  (39)  and  (40)  are  now 


e  «  fuu  ♦  1/2  fuuu< 


u'2  ■  1/2  a  u2  ♦  1/6  g  u3 
'  *uu  6uuu 


The  coefficients  f  and  gyu  are  identical  to  those  in  Reference  1,  and  the 

other  two  coefficients  are  only  slightly  modified  so  the  behavior  of 
solutions  is  still  the  same  as  described  in  Reference  1.  Thus  if  c  is  chosen 
such  that  g  >0,  solutions  are  solitary  waves. 


Ill 


Clearly  e  must  be  chosen  such  that  |3«uu|  <|  guuJ  •  Propagating  bulges 
(necks)  occur  if  guuu  <  0  (suuu  >  0).  The  situation  for  a  propagating  bulge 
is  shown  in  Figure  6. 

In  terms  of  the  elastic  moduli  of  linear  elasticity,  the  second 
derivatives  of  W  used  above  may  be  written  (see  Reference  1) 


Wff  *  X  +  2y,  Wuu  -  4(X  ♦  y),  Wuf  -  2X,  Wg  -  1/4  a2y  (52) 

where  X  and  y  are  the  Lame  elastic  constants.  Therefore 


.  16  (1  ♦  y)  (E  -  pc2) 

u  a2(X  +  2y  -  pc2)  (y  -  pc2) 


where  E  is  Young's  modulus.  The  conditions  on  c  will  be  satisfied  if 
0  <  ((  pc2/E)  -  1)  «  1. 

If  g(0)  =  u'  2  j*  0,  then  (50)  must  be  replaced  by 

u'2  B  "a*  *  V  +  8uu“2  +  1/6  ®uuuu^  (54) 

2 

In  this  equation  it  turns  out  that  if  jau^|  ■  6«1,  then  gu  =  0(6  ),  and  both  gyu 

and  guuu  differ  from  their  previous  values  by  terms  that  are  0(62).  If  c  is 

chosen  so  that  g  is  negative,  then  solutions  to  (54)  are  periodic.  The 

situation  that  is  envisioned  here  is  shown  in  Figure  7.  Solutions  will 
oscillate  between  the  points  u^  and  u g,  and  may  be  expressed  in  terms  of 

elliptic  integrals.  (See  Reference  6,  Chap.  17  for  the  appropriate  transfor¬ 
mations.)  Because  of  the  restriction  that  g^u  <  0,  not  all  speeds  are 

possible  for  periodic  waves.  In  the  linear  limit  this  gives  the  familiar 

2  2 
conditions  y  <  pc  <  E  or  X  ♦  2y  <  pc  . 

VI.  SUMMARY  AND  DISCUSSION 

It  has  been  shown  that  there  is  a  complete  mathematical  analogy  between 
theories  for  a  nonlinear  rod  and  a  granular  material.  Since  rods  are  familiar 
objects  to  applied  mechanicians,  but  study  of  granular  materials  is  of  rela¬ 
tively  recent  origin,  the  analogy  should  be  useful.  For  example,  it  was  noted 


6M.  Abramauitz  and  I.  Stegun ,  eda. ,  Handbook  of  Mathematioal  Functions 
Dover  (1965). 


that  a  characteristic  length  of  /2k  could  be  identified  for  a  granular  mate¬ 
rial,  but  even  more  important  any  statement  about  the  dynamics  or  deformation 
of  a  rod  has  an  exact  counterpart  for  a  granular  material.  Thus,  for  static 
homogeneous  deformations  of  a  granular  material  the  equation  of  balance  for 
micromomentum  (2)  reduces  to  g(v,o,e)  *  0,  which  defines  an  equilibrium 
relation  v  =  v(e),  and  the  negative  of  the  slope  of  this  relation  defines  the 
analog  of  Poisson's  ratio.  Furthermore,  along  this  equilibrium  line  there 
will  be  a  modified  stress-strain  relation,  o  =  o  (e),  the  slope  of  which 
corresponds  to  Young's  modulus  E.  In  fact,  it  is  Young's  modulus  that  appears 
as  the  first  term  in  the  numerator  of  (47).  For  the  granular  material  the 
analog  modulus  may  be  written 


e  e  -  e 
vv  ee  ve 


Since  it  is  known  that  Young's  modulus  has  fundamental  importance  for  rod 
deformations,  the  analog  modulus  must  have  fundamental  importance  for 
deformations  in  granular  materials.  In  particular,  for  low  frequency 
disturbances  with  wavelengths  substantially  longer  than  the  characteristic 
length,  it  is  to  be  expected  that  the  analog  modulus  will  dominate. 

It  has  been  pointed  out  that  both  rods  and  granular  materials  support  two 
kinds  of  shock  wave.  One  is  the  familiar  longitudinal  shock  for  both  cases, 
but  the  second  is  somewhat  unfamiliar.  The  second  type  of  shock  is  associated 
with  the  micromomentum  and  is  accompanied  by  a  discontinuity  in  the  gradient 
of  radial  strain  for  the  rod  or  by  a  discontinuity  in  the  gradient  of  volume 
fraction  for  the  granular  material.  In  the  rod  such  a  discontinuity  is 
associated  with  a  radially  symmetric  shearing  motion  so  it  can  be  readily 
visualized,  and  it  is  not  surprising  that  the  speed  is  controlled  by  the  bulk 

p 

shear  modulus,  w  /%a  ,  which  is  just  the  second  Lame  constant  u  for  a  linear 
4Q 

material.  The  analog  for  the  granular  material,  however,  is  totally  unfamil¬ 
iar,  but  may  be  of  some  importance.  The  fact  that  the  stored  energy  is  even 
in  q  (or  v' )  has  also  been  pointed  out.  The  symmetry  of  the  amplitude  curves 

for  q+  :  0  in  Figures  2-5  and  the  slight  biasing  to  one  side  or  the 
other  for  q*  a  0  are  direot  consequences  of  that  fact. 

Finally,  it  has  been  pointed  out  that  finite  amplitude  periodic  waves  and 
solitary  waves  can  propagate  in  either  a  rod  or  a  granular  material.  There 
are  two  spectral  branches  for  wave  speeds,  as  was  noted  at  the  end  of  the  last 
section.  Detailed  calculations  with  the  elliptic  integrals  that  satisfy  (54) 
would  show  the  effect  of  finite  amplitude  on  wave  speed  for  periodic  wave- 
trains  with  a  fixed  wavelength. 
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